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ABSTRACT 


In a time-division multiple-access (TDMA) system, many 
transmitting stations time share the same RF spectrum. 
Individual transmissions must be identifiable to distant 
receiving stations and must be synchronized so that no 
interference overlap occurs. A unique word sequence can 
serve both purposes. | 

The basic concepts of a satellite communications TDMA 


system are briefly reviewed, with emphasis on functional 


receiver reguirements. A correlation detection scheme is 
hypothesized, from which false and miss-detection 
probability expressions are derived. Criteria are then 


established for the selection of sequences suitable for use 
as unique words. Representative unique words of length 15 
through 20 are included, along with an explanation of the 
difficulties encountered in devising an efficient comput 
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I. INTRODUCTION 


The possible use of time-division multiple-access (TDMA) 
in satellite communications has received a considerable 
amount of attention jn recent years. In a TDMA system, many 
stations use the same satellite repeater on a time shared 
basis. Transmitter timing must be synchronized so that only 
one station's signal accesses the repeater at any given 
time; and, further, so that successive Signal bursts are 
closely spaced in order to make efficient use of the 
repeater capacity. Continuously changing ranges from ground 
Stations to the satellite (Appendix B) make transmitter 
synchronization a non-trivial problem. Once the 
Synchronization problem is solved, TDMA offers several 
advantages over current  frequency-division  multiple-access 
procedures 11,2% Schmidt [3] outlined the proposed 
application cf TDMA to the INTELSAI-4 series commercial 


communication satellites. 


À. TDMA SYSTEM DESCRIPTION 

As each station must transmit in timed bursts, 
transmissicns are necessarily digital and provisions must be 
made for the sampling (and multiplexing) of the transmitted 
information channel(s). This process must be timed so that 
it is in consonance with the burst repetition frequency 
(BRF) and transmitted data rate. When analog-to-digital 
conversion 1S involved, the analog sampling rate must be 
equal to, or a multiple of, the BRF; a sample storage 
capability is rejuired when the sampling rate is a multiple 
of the BRF. When information to be transmitted is in the 
form of a continuous bit stream, a buffer is required to 
Store a given portion of the bit stream for transmission at 
the (higher) data transmission rate. In either case, the 
number of information channels that can be transmitted in a 


burst of given length is a direct function of the ratio of 





transmission rate to information rate. Time-division 
multiple-access, then, is in effect a higher level fcrm of 
time-division multiplexing. TDMA is a time sharing at RF, 
while TDM is a tjme sharing at baseband, before carrier 
modulation. The process of several transmitted signals time 
Sharing a satellite repeater is conceptually similar tc that 
of several information channels time sharing the same RF 
carrier. 

Given that transmitted Signals are neatly interleaved in 
their arrival at the satellite, they will be retransmitted 
in the same fashion and received on the ground without any 
interference overlap. A receiver must first recognize when 
a burst is of interest and then, in the case of multiflexed 
transmissions, select particular information channels within 
the burst and route them to appropriate end terminals. In 
order to do so, the receiver must be provided with certain 
control information. 

1. Receiver Functional Requirements 

SuccesSive received bursts have been originated by 
Gifferent, in general Widely separated, stations whose 
carrier frequencies are not coherent in frequency and phase. 
Therefore, prior to the demodulation of each burst, the 
receiver local oscillator must be locked in frequency and 
phase to that of the incoming carrier using phase-locked 
loop techniques. (Use of differential phase-shift keying 
removes the carrier recovery requirement at ihe expense of 
error-rate performance (8S ].) Similarly, once carrier 
recovery has been accomplished, the receiver bit tining 
oscillator (clock) must be synchronized to the incoming bit 
stream. Bitzer L5] proposed a system design for 
establishing coherent bit timing among all Stations and 
Suguri, Doi and Metzger [4] reported on the successful test 
of this type of system. In such a system, receiver and 
transmitter clocks are continuously aàjusted to maintain 
coherency within acceptable limits, obviating the need for 
gross adjustrents to the receiver clock preparatory to burst 


demodulation. 









With the local oscillator and clock synchronized, 
demodulaticn of the incoming signal occurs. Enbedded 
somewhere in the beginning of the burst must be information 
which the receiver can extract to obtain transmitting 
station identity, word (channel) synchronization and channel 
status. Demodulation/demultiplexing of the information 
channels ccmmences thereafter. 

2. Irass and Durst Format 

Satellite repeater access time is divided into 
frames, a frame being the length of tine allotted for the 
transmission of one burst from every station. The frame 
length is typically 125 usec, or a multiple thereof, to 
conform tc the Nyquist Sampling rate for voice (8KHZ2) [3,6]. 
A station is assigned a time slot within the frame for the 
transmission of its burst, the length of which depends on 
the anticipated channel requirements. A guard time of 
typically 200 nsec or less separates the time slcts to 
preclude burst overlap. ; 

Each burst commences with a preamble, follored by 
the information channels. A typical format is depicted in 
Fig. 1-1. <A fixed pattern of dummy bits is inserted at the 
beginning for the purpose of locking up the local oscillator 
and clock. A constant string such as 00009eee for carrier 
recovery, followed by an aiternating pattern such as 1010ese 
for bit timing recovery is typical. Following the recovery 
bits is a fixed sequence of bits unigue to the transmitting 
station, the "unique word." More vill Le said about the 
unique word below. The remainder cf the preamble contains 
control signals and orderwire transmissions necessary for 


assignment and control of the information channels. 
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Figure 1-1 
Frame an Burst rormat 


TDMA efficiency is a measure of the total tine 
allotted to the transmission of actual inforuation versus 
the total available repeater time. The preamble is wasted 
transmissicn time in this sense and considerable effort has 
been expended to reduce its length [4,5,8,9]. 

3. Unigue Word Functions 

To the extent that it is recognizable tc tne 
receiver, the unique word can be used to identify the 
Originating station. It may also be used as a timing 
reference to establish word synchronization, resetting a bit 
counter which properly routes the control and orderwire 
portions cf the preamble and indicates the beginning cf the 
first information channel. 

It has been assumed up to this point that ail 
Fransmittang stations are synchronized so that no burst 
overlap SC UTS. Nothing has been said about how 
synchronization is acquired or maintained. The usual method 


of doing sc is to use the unique words as timing marks 
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within the frame. One station is designated as the master 
station and its unique word used as the fundamental timing 
reference for the frame. All other stations compare the 
tine difference between their own unique word and the 
reference, as recejved from the satellite, and adjust 
transmitter timing as necessary, one bit per frase, to 
maintain proper separation. This problem of transmitter 
Synchronization is of interest herein only in that the 
unique word plays a vital role, and will not be further 
pursued. References Zeck and 10 offer a detailed 
description cf the various factors which must be considered 
and provide insight as to how a synchronization package 


might be instrumented. 


B. DETECTION AND SELECTION OF UNIQUE WORDS 
A TDMA system of the type outlined above can be 
operated successfully only if the receiver is able to detect 
the receipt of a unique word and differentiate one from the 
other. Some type of correlation detector is infcrred. 
Proposed detector design features and operation are 
presented in Section II. Then, in Section III, criteria are 
established for the selection of unique words which yield 
acceptable miss<detection and false detection probabilities. 
Interest in the subject was motivated primarily by the 
discussicns of Schrempp and Sekimoto [11]. Much of the 
detector design philosophy is borrowed from their work, and 
their Approach to unique word selection has been 
instrumental in obtaining an appreciation for the 
fundamental requirements which must be considered. The 
selection criteria set forth in Section III and the 
subsequent computer search results represent the major 


products cf this research effort. 
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II. DETECTION OF UNIQUE WORDS 


AME 2 ap eg uU s EP - u 


A unique word is a fixed binary sequence of length N. 
Given a knowledge of the elements of the sequence, detection 
can ke accomplished by comparing received bits, N at a time, 
with a replica of the sequence stored in the receiver; i.e., 
by using correlation techniques. The probability of missing 
the unique werd (miss~detection) and the probability of 
receiving the sequence at random (false detection) are 
fundamental considerations in the evaluation of detector 


effectiveness. 


A. BASIC DETECTOR DESIGN 

A digital correlation detector of the type depicted in 
Fig. 2-1 is postulated. The element values of the unique 
word are set senj-permanently in a register and the 
demodulated kit stream is clocked through the shift 
register. The corresponding elements of the two registers 
are added, mcdulo-2, once each bit period and the results 
are summed algebraically for input to the threshold network. 


The ncdulo-2 operation is defined by, 


060=0 
00 1-7 1 
160-7 1 
161-0 


And tbe input to the threshold network is, 


H 
nou» Se DEUS) (25%) 
i-1 2 1 
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R therefore equals the number of disagreements between S 
and  U. The threshold network is designed so that a 
detection pulse is produced whenever there are h or fewer 
disagreements; that is, whenever RSh. The setting of a 
threshold permits the detection of unique words which 
contain up to h transmission errors and thereby reduces the 
niss-detection probability. The probability of false 
detection increases in doing so, however, due to the 
increased opportunities for a random match-up. 

It shculd be noted that R, as mathmatically defined, is 
in fact the inverse correlation function of S and U. When 
the elements of S and U all match, R-0; and when there is 
complete disagreement,  R-N. Alternative detector designs 
can be envisioned which would produce the true correlation 
BEunctron at R; e.g., construct adders that produce the 
complement of the modulo-2 sum and design the thresnold 
network so that a pulse is produced whenever R2N-h. Of 
course, in the end, the voltage levels associated with  '0' 


and 171 are wnat determine the actual "value" GE Cen 








Whichever of the various alternatives are chosen, the end 
result is the same: a detection pulse is produced whenever 
there are h or fewer disagreements between S and U. 

The capability of receiving more than one transmitting 
station is provided by a parallel bank of detectors, each 
set to the unique word of one of the transmitting stations. 
The received bit stream 1s passed Simultaneously through 
each detector and the resulting detection pulses are used to 


control receiver functioning, as described in Section I.A.3. 


B. STEADY-STATE APERTURE GATING 

The probability of false detection is reduced 
Significantly if the detector can be gated so that it kegins 
looking for a unigue word just prior to its receipt, and is 
otherwise inactive. The detection pulse itself may be used 
to generate such a gating signal under steady-state 
conditions; that is, once the unique word has been initially 
acquired. 

A unique word is iransuitited every T seconds, once each 
frame, except when a correction is required to maintain 
proper position within the frame. If corrections are made 
at the rate of one bit per frame {9,10}, the elapsed time 
between unigue word transmissions is bounded by T+b, where b 
is the period of one bit. The time separation of successive 
received unicue words is further affected by relative notion 
between the satellite and the respective transmitting and 
receiving stations. Let d, in bits per second, te tle 
maximum expected doppler induced drift of successive unique 
word detection pulses. The maximum doppler shift in one 
frame period is then 4Tbd and the elapsed time between 


Successive receipts of a unigue word is bounded by, 


ote Gb» fad OG) 


T 
I 


H 
tt 


T + b(1 + Tà) (2.2) 
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This relationship can be expressed in terms of bit 
periods ky dividing through by b. Doing so establishes a 
limit on the number of bit periods separating successive 
receipts of a unigue word. 


N 
r 


T /b = T/b X (1 + Td) 
E 


N 
r 


n + (1 + Td) 


where n is the number of bits in one frame. 
Defining D as the maximum displacement (in bits) c£ the 


unigue word from its "periodic" position, 
D = + Td (2.3) 


N =n tr D (2.4) 


Upon receipt of a unique word, the detector is turned 
off, gated back on T-bD seconds later, gated off by the next 
unique word, etc. raise detections can therefore result 
only from a random match-up within the aperture immediately 
preceding the unigue word. 

The guestion remains as to what to do when the unique 
word is missed due to errors occurring in transmission. 
Once the detector is gated on, is it left on until a 
detection pulse is received (in which case a missed unique 
word causes an immediate transistion to the acquisition 
mode); or, is the detector gated off at time T+bD and then 
gated pack on at time 2T-bD, one frame later? In the 
latter case, how many successive misses are required kLefore 
cancelling the aperture and reverting to an acquisition 


node? The answer depends for the most part on particular 


It is assumed that the aperture width remains constant 
and would not be expanded as successive misses occur. 
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system specifications. When a unique word is missed, all 
informaticn contained in the burst is lost. Having missed 
it once, that is having received the unique word with more 
than h errors, there is a given probability that it will 
contain more than h errors the next time. Calculation of 
the probability of a unique word containing more than h 
errors two or more times in succession is not enough 
however. The doppler induced burst drift rate must also be 
considered. 

If the drift rate is such that the unique word is 
received outside the aperture, it will be missed, even if 
error free. On the ether hand, reversion to an acquisition 
mode may involve a considerable penalty in lost information 
due to a long average acquisition time. (See Section II.C 
below.) Thus, a desire to avoid the acquisition mode must 
be balanced against the diminishing probability (as 
Successive misses occur) of receiving the unique word in the 
aperture mode. AS stated earlier, this is a systens 
dependent decision, and one that can best be reached 
empirically in an operational environment. For present 
purposes, the rather safe assumption will be made that at 
least two successive misses are required before a transition 
from the aperture to acquisition mode is effected. À 
detection pulse turns off the detector and sets a  ccunter 
Which gates the detector on T-bD seconds later. If no 
detection pulse occurs, the counter turns off the detector 
at time TtbD. The aperture is of width 2D+1 (bits). 

The maximum burst drift which may be expected when using 
a geosynchrcnous satellite is computed in Appendix B. 
Equation B.5 can be substituted into equation 2.3 to obtain 
values for the maximum displacement of a unique word over 
one frame period. Table T2=] contains representative 
calculated results for the shipboard platform discussed in 
Appendix Bas well as for a platform with a relative 


velocity on the order of 1000 knots. 
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v(nzsec)  f[(Mbps)  d(bits/sec)  Tí(nsec) D(bits/franme) 





-5 

61 1 0.407 0.125 1+5.1x10 
61 1 0.407 3.0 1.001 

6 1 50 2043 0.125 1. 003 

61 50 20.3 3.0 1. 06 

61 500 203 0. 125 1.03 

61 500 203 3.0 1.61 

2000 1 13.3 0.125 1.002 

2000 1 13.3 3.0 1.04 

Table 2-1 


Bit Rate Doppler (d) and Haximum Displacement (D) 
of a Unique Word Over One Frame Period for Given 
Values or Stataon/Satellite Relative Motion (Vv), 
Bit Fregquency(f) and Frame Period (T) 


C. ACQUISITION GATING 

The steady-state operation of a unigue word detectcr was 
described in the preceding section. Here the process of 
acguisition is considered. When a receiver is activated, 
each detector must acquire the associated unique word within 
the TDMA frane. Likewise, a detector must be switched to 
the acquisition mode whenever steady-state synchronization 
is lcst. Acquisition, or re-acquisition, may be effected by 
sinply passing the received bit stream through the detector 
until a detection pulse is received. Then, as in the 
steady-state mode, the detection pulse switches the detector 
off and starts an aperture gate counter which in turn 
Switches the detector back on in the next frame. If a 
second detection pulse is received within the aperture, the 
detector is transferred to the steady-state mode. 
Otherwise, the acquisition process continues. 

When in the acquisition mode, the bursts of all active 
stations are being passed ccntinuously through the detector 


and the probability of false detection is thus relatively 
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nag hi. Should a false detection occur, a false aperture is 
set and the next actual unique word detection pulse will be 
missed. Even when the detection is valid, there is a 
non-zero probability that the unigue word will be missed in 
the next frame. In either case, acquisition, and transfer 
to the steady-state mode, is delayed for at least two frame 
periods. The actual acquisition time may thus span several 
frame periods. 

It is possible to construct a detection pulse storage 
network which would look at detection pulse locations over 
two or more frame periods and use majority logic to set an 
aperture corresponding to the true detection pulse. 
Although acquisition time would be improved, implementation 
could prcve relatively costly because each detector would 
require its own network. An alternative is proposed, which 
makes use of the same aperture network employed in the 
steady-state mode to gate the detector and reduce the 
probability of false detection. 

1. Frame Reference Gating 

Each burst has an assigned position within the TDMA 
frame. It is known that a particular burst will be received 
a given number of bits after the receipt of the reference 
unique work, plus or minus the guard band cushion. ` The 
guard band is set so that it prevents burst overlap; that 
is, it represents the maxinum expected displacement of a 
burst (or unigue word) from its nominal position. It is 
therefore proposed that, during acguisition, the 
apertureycounter network of a detector be keyed by the 
receipt of the reference unigue word. The aperture is set 
so that the detector is gated on at the assigned position of 
the corresponding unique word, minus the number of bits in 
the guard band, B. If the unique word is missed, the 
detector is switched off B bits after the noninal unique 
word position. (The aperture width is 2B*1.) If the unique 
word is detected, the aperture is set for the next frame. 
Two consecutive detections Key a transfer to the 


steady-state mode. 
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In a demand access system, where at least some of the 
Stations are not assigned a permanent time Slot, use of 
frame refererce gating is feasible only if all stations can 
be provided with time slot activation and deactivation 
informaticn, or can otherwise predict where a given 
station's burst may appear within the frame. Whether cr not 
it is possible to provide advance activaticn notification to 
all stations depends on how a part time station gains access 
to a time slot. If access is centrally controlled 
(automatically), then assignmsont signals from the 
controlling station can be copied by all stations and used 
to set the acquisition gate. If a station gains access by 
seizing a vacant time slot, a receiving station must have 
other means of predicting where the burst will occur. The 
results of the computer search for unigue words, reported in 
Section IV; suggest a procedure for meeting this 
reguirement. 

For a given unique word length, several different sets 
are available which have good cross-correlation properties. 
Unique words from any one set could be assigned to certain 
Stations shich collectively have exclusive use of a block of 
consecutive time slots within the frame. As many such 
blocks as necessary could be made available to different 
groups cf stations.. The number of time slots assigned to 
each group wculd be based on collective useage requirements, 
and could be varied on a scheduled or demand basis. (Demand 
access implies that the number of time slots would be less 
than the number of stations.) Those receiver detectors 
corresponding to the demand access transmitters could then, 
when in the acquisition mode, be gated on at the beginning 
of each burst within the demand access  sub-frame. The 
probability of false detection will be somewhat higher than 
for the single gate procedure proposed above; but, with low 
cross-correlation between unique words of a set, will be 


much lower than for "wide open" acquisition. 


21 





Having cutlined the basic detector design concefts, a 
basis has been formed for the selection of sets of unique 
words. A detailed analysis of the attendant detection 
prokabilities is presented in Section C. Simply stated, the 
goal is to find unique words which yield low niss-detection 
and false detection probabilities. Criteria will be 


established in this section for doing so. 


A. LENGTH 

The fewer bits in a unique word, the lower the 
miss-detection probability. Conversely, the longer the 
unigue word, the lower the false detection probability. The 
use of a detector threshold and aperture nullifies to sone 
extent the importance of these opposing considerations, but 
in a gross sense they can never be eliminated as factors 
influencing the ultimate choice of unique words. In fact, 
length and the threshold setting are the cnly parameters 
available for altering the miss-detection probability. 

From the standpoint of efficiency, it is desirable to 
keep unigue words as short as possible. The longer each 
station’s unique word, the less time is available for the 


transmission of information. 


EXA UI TOCORRELATION CRITERIA 

Operation of the detector in the steady<state mcde is 
described in Section II.B. Loss of channel inforaation will 
result if (1) the unique word is missed, or (2) a false 
detection occurs in the aperture prior to the uniyue word. 
Either may lead eventually to a drop-sync condition. As (1) 
and (2) are equally damaging, it would seem appropriate to 
select unique words in such a tay that their combined effect 
is minimized. Under the conditions described in Section 


II.B, it is largely a matter of unguantifiable chance 
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whether a miss-detection or a false detection (within the 
aperture) is the least desirable occurrence. Accordingly, 
the following condition was established: if the received 
unique word contains h or fewer errors, the false detection 
 probability must be zero within the aperture. A false 
detection pulse will not be produced unless the unique word 
is going to he missed anyway. 

As it is possible for a unique word detection pulse to 
occur on the last cycle within the aperture, and as the 
aperture is of width 2D*+1, the foregoing condition can be 
met only if the probability of false detection is zero for 
the 2D cycles preceding the unigue word. This condition can 
in turn be satisfied only if the autocorrelation function of 


the unique werd is such that, 
R(t) « N- 2h - t ; 1<t<2D (3. 1) 


where t is an integral time variable, the basic unit of 
which is cne bit perjod, and, 


N-t 
R(t) wt U (i)6U (i*t) 
1 


This condition is presented graphically in Fig. 3-1. 
R(t) is egual to the number of agreements hetween the 
(error-free) N-t displaced unique word bits and the detector 
setting. Up to h errors are permitted in the unique word 
under the prescribed conditions and it is possible that 
these errors may be such as to produce additional 
agreements, giving a total of R(t)+h agreements frcm the 
displaced unigue word alone. There are t non-unigue word 
bits in the detector. If all of them match with the 
detector setting, a combined total of R(t)+h+t agreements 
results. As N-h or more agreements are required to produce 
a false detection, the probability of such occurrence is 
zero if R(t)th+t<N-h, and the validity of inequality 3.1 is 
established. 
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In the following section, a condition for the 
cross-correlation function will be established. This same 
condition must be met by the autocorrelation function 
outside the aperture in order to meet specified limits on 
the probability of false detection during acquisition. It 
is listed here, in advance, in order to complete the 
autocorrelation specjfications. 

R{tje< N- h= t ; 2D<t<N—h (3. 1a) 
R(t) = 0 ¢ N~hsSt<N~1 (3.1b) 

There are 2" sequences of length N from which to choose 
unique words. However, symmetry considerations can be 
invcked tc reduce the number which must actually be checked. 
First, because the mod-2 summing operation is commutative, a 
sequence and itS complement, obtained by inverting the value 
of each bit, have the Same autocorrelation function. 
Second, because autocorrelation is an even functicn, a 
Sequence and its transpose (or mirror), obtained by 
inverting the order of the bits in the seguence, have the 
same autocorrelation function. Therefore, once a sequence 
is checked it iS unnecessary to check its complement, 
transpose or complement of its transpose. The numrer of 
required checks is not reduced by a factor of four, as it 
would appear at first glance, however, because the transpose 
and corpiement transpose pair may be identical tc the 
original sequence and itS complement. The required number 
of checks can be determined by consideration of the effect 
of the complement and transpose operations. Although the 
results were not used directly in the computer search 
routine, they do provide an indication of search efficiency, 
and are presented here as a matter of peripheral interest. 

In crder for the transpose or complement transpcse to be 
identical to the original sequence, the second [N/2] bits of 
the sequence must be the transpose or corplement transpose 
of the first [N/2] bits, vhere [N/2] equals the largest 
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integer less than or equal to N/2. (If N is odd, the middle 
bit of the sequence does not change position under the 
transpose operation.) The first [N/2] bits may be arranged 


BZ, 


in different ways, each arrangement having an 


associated transpose and complement transpose in the second 


Se? 


N 
half. Therefore, 2x of the original 2 sequences have 
redundant (identical to original sequence) transposes or 
complement transposes, and, for this subset, the required 


nunber of checks can be reduced by a factor of two only (by 


pairing a sequence with its complement). There are 


N [N/2]*1 Wei 
22 seguences remaining to be considered. These can 


be grouped by fours (sequence, complement, transpose, 
complement transpose) to reduce the number of reguired 
checks by four. The total number of sequences of length N 
which must Le checked in search of unique words meeting the 
autocorreiaticn coudiiion is tüerefore, 


N N/2]*1 
A ` ëch 


N/2 
Required Checks + al AE 


N=2 N/2 3-1 
2 + al ed 


The computer program used to accomplish the 


autocorrelation sort is listed in Appendix D. 


See ROSS-~CORRELATION CRITERIA 

The cross-correlation between the unique words used by 
stations sharing the TDMA frame is a crucial factor 
affecting the false detectic probability during (ungated) 
acquisiticn. A cross-correlation characteristic that yields 
a high probability of false detection is particularly 


unsatisfactory because of the periodicity of its occurrence. 
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In such cases, the detection aperture may be falsely set 
time and time again over several consecutive frames. 

The probability expressions derived in Appendix C were 
used as a basis for the establishment of the following 
cross-correlation condition: the probability of false 
detection must be zero unless the received unique word 
contains cne or more errors. This condition also applies to 
the head and tail of the autccorrelation function, outside 
the aperture. The cross-correlation condition is thus 
similar to the autocorrelation condition, differing only in 
that the unigue word must be error free, vice h or fewer 
errors, to guarantee zero false detecticn probability. 
Using R(t) again to represent the correlation function, the 
conđition cn R(t) is (analogous to inequality (3.1)) 
R(t)<N-h-Itt. AS the correlation function is always 
non-negative, this condition can be met only when [t{<N-h. 
For higher values of |t| the condition must be relaxed to 
R(t)=0. Summarizing, sets of possible unique words, having 
passed the autocorrelation test, are retained oniy if the 


cross-correlation between all words of the set satisfy, 
R(t) <N- h - Iti ; [t{<N-h (3. Za) 


EI) (0 :  N-hs[t|XN-1  (3.2b) 

The ccnditions placed on R(t) can be used to further 
refine the probability expressions of Appendix C. It will 
be assumed that bits of fixed value (at least N-1 of them) 
precede the unique word and that bits of random value 
follow. This assumption is based on the system 
specifications outlined in Section I, and, although not an 
absclute necessity at this point, is made in order to reduce 
the number of situations which must be considered. Other 
Situations, when and if occurring, may he handled in a 
completely analogous fashion. 

First, the probability of false detection in the fixed 
bit regicn preceding a unique word (a unique word other than 


the one for which the detector is set) will ke considered. 
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The point where the unique word is completely contained in 
the detector defines t=0. At t=- (N-1) the first bit cf the 
unique wcrd enters the detector.  Equation(s) C.8 apply for 
~ (N-1}) st<0, with m=N. As in Appendix C "number of 
agreements" refers to matches between the detector setting 
and the value of the fixed bits as transmitted. The number 
of agreements, k in equation C.8, is equal to R(t), the 
number of matches within the unique word bits, plus the 
number of matches within the remaining |t| bits. If all of 
the latter match, 


k = R(t) ¢ jt] 


Therefore, the cross-correlation condition 3.2a 


(R(t) <N-h-[{t]) 1S equivalent to the equation C.8 condition, 
k * h « N 


The cross-correlation condition 3.2b in itself does not 
quarantee that k*h«N. It is necessary to impose the 
additional systems dependent condition that fever than  N-h 
agreements can be permitted between the non-unigue word bits 
for - (N-1)<t<- (N-h). Doing so, the probability of false 


detection at any given step is (from eguation C.8e), 


Pí(^(t)] - 


N-h-k h*k 
C(N-k,N-h-k) p (1-p) >; -(N-1)sts0 (3.3) 


if k takes on its maximum value (N-h), the probakility 


of false detection becomes, 


N-1 
P{A(t)} = (h+1)p(1-p) "e (N-—1)Sts0 (3.4) 
max 


The probability of false detection in the randcn bit 
region following the unique word is also given by 


equation(s) C.8. In this case m equals the number of unigue 
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word bits in the detector (there are no other fixed bits) 


and the number of fixed matches, k, equals R(t). 


w = N - ft] 
REESEN 
Therefore, the cross-correlation condition 3.2a 


(R(t) <N-h-|]t]) is equivalent to the equation C.8 condition, 
k + h < on 


And the probability of false detection at any given step is 


given by equation C.8c (with m and k given by above). 


P{A(t)} = 


= -k- 


m m h kth 
C (n-k,m-k-h) p (1- p) ; OQO<t<N<h (3:5) 


N 
(1/2) 
When k takes on its maximum value (R(t) =N-h-t-1), 


equation 3.5 becones, 


o Y Am 


t N-t-1 
P {A (t)} = (1/2) (h*i)p(1-p) ;  O«t«N-h (3.6) 
max 


For N-hsts<N-1, m=h,h-1,h-2,00e,1 and k=R(t)=0, giving 


m = N- t 
and, 
ont eee! 
The probability of false detection is (from equation 
ED»), 
i= m _hrn 
EU (v3 — (ph ën Eifel (3.7) 


t N-t t 
(1/2) (7p) — 2, — C(t,J)i EA 


30 





Equation 3.7, unlike equations 3.3 and 3.9 is 
independent of parameter variation [i.e. independent of k, 
R(t) j] and thus "maximum P{A(t)}" has no meaning. 

Equations 3.3, 3.5 and 3.7 may be used to calculate the 
probability of false detection as each unique word in the 
frame is passing through the detector. Equations 3.4, 3.6 
and 3.7 frovide a worst case limit for these probabilities 
(when unique words are chosen using the cross-correlation 
criteria 3.2a and 3.2b). These individual fixed-bit event 
probabilities are the ones needed, among others, in equation 
C.10 to oktain an upper bound for the probability of one or 
more false detections (F) over the period of one frame. Or, 
over a narrower time span, they can be used to compute an 
upper bcund on the probability of one or more false 
detections during the 2N-1 detection events surrounding a 


given unique word; that is, to compute, 


F! provides a means of obtaining a relative measure of 
the goodness of different unique word choices without having 
to resort to the more lengthy calculations required by 
equation C.10. Substituting from equations 3.4, 3.6 and 
B7, 


EN - 


N N-h t N-t-1 
F' € N[ (h¢1) p(i-p) ] E (17/2) (h*!)p(l-p)j 


aem 


Bst. ee 


Ec (1/2) (1-p) Aë c (t, 3) (3.8) 

The computer program used to check cross-correlaticn is 
listed in Appendix D. Each sequence which passes the 
autocorrelation test 1S cross-correlated with every other 
such sequence, yielding pairs of acceptable unique words. 
In crder to form maximal length sets of acceptable unique 


words, it is necessary to perform an additional matching 
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routine. For example, if UN 1 cross-correlates well 
(matches) with UW2, UW3 and UW4&, then a set of three is 
obtained if UW2 also matches with UW3 or DND, or UW3 matches 
with UW4; a set of four is obtained if UW2 matches with both 
UH3 and UW4 and UW3 matches with JUW4. Various matching 
algorithms were programmed, none of which performed ina 
completely acceptable fashion. At N=17, the number of 
checks which must be performed becomes so great that, for 
the algorithms used, computation time runs into days. In 
order to oktain results for N=17 and higher, various 
expedients were employed, as described in the following 


section, to reduce the required computer time. 
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IV. PRESENTATION OF RESULTS 


The SDS-9300 digital computer was used to perfcrm a 
correlation sorting routine patterned after the criteria 
established in the preceding section. At the outset, it vas 
realized that a significant amount of machine time would be 
required to cbtain results; but, even so, the magnitude of 
the problem Was not fully appreciated until program 
execution was initiated. The SDS-9300 suas chosen over the 
IBM-360 kecause of the local availability of multi-hour 
blocks cf computer time, realizing that its limited on-line 
memory (32K) might prove restrictive. 

The choice of detector parameters and the essential 
findings of the computer search effort are covered in the 
following paragraphs. Selected sets of unigue words are 


listed in Appendix A. 


A. CHOICE OF DETECTOR APERTURE AND THRESHOLD 

The  rotivation for use of a steady-state aperture 
technique is provided in Section II.B. As the aperture 
width must be at least 2D+1, where D is the maximum 
displacement of a unique word in one frame period, Table 2-7 
leads one to a choice of a four or five bit width. An odd 
width is referred in order that the aperture counter can be 
run directly off the receiver clock; i.e., so that detector 
activaticn/deactivation occurs at the beginning rather than 
in the middle of a bit cycle. An aperture width of five 
bits is thus the mimimum acceptable, and yet at the sane 
time provides sufficient margin to cover unique word drift 
over more than one frame in the event of a miss-detection. 

The choice of a detector threshold (h) 1s determined Dy 
the setting cf an acceptable limit for the miss-detection 
Peenacitity (0), Which in turn is a function of the bit 


error rate (I). For this purpose, p uas assigned a value of 
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10 [11] as being the worst-case bit error rate that could 


be tolerated on a satellite communication channel. 


Equations C.1, C.2 and C.3 were used to compute Q and the 
mean time between miss-detections (A) for various values of 
N and h. The results are listed in Table 4-1. The tines 
given for HM assume a frame period of 125 usec. 

Acceptable values for li occur for h-2 and higher. The 


fact that the probability of false detection increases with 
h dictates the selection of the minimum acceptable 
threshold; viz. h-2. 


NO h Q n 
i5 0 EES 83 msec 
1 1.05x10 ^ 2 min 
2 Mo 4 76 hrs 
3 1.36x10 "7 29 yrs 
17 0 1.70x10 ^ 74 msec 
A conos. 92 sec 
2 6.79x10 | 51 hrs 
3 2.38x10 - 17 yrs 
20 0 2.00x10 ^ 62 msec 
4 1.90x10 ^ 66 sec 
2 1. 14x102, 30 hrs 
3 4.84x10 8 yrs 
Tanie 4-1 


Miss-Detection Probabilities 
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B. UNIQUE WORD SYMMETRIES 

Certain characteristics of the sequences which pass the 
autocorrelation test were used to simplify the 
cross-correlation sort. These same characteristics also 
provide a degree of flexibility in applying the results to 
particular operating systems. 

A seguence can pass the autocorrelation test only if its 
autocorrelation is zero for the first and last three values 
at the head and tail of the distribution (with h=2). 
Therefore, the first three bits of the sequence must te ali 
zercs or all ones, and the last three bits must all te of 
Opposite polarity. The need to check less than one-half of 
the total possible sequences was explained in Section III.B. 
AS the autocorrelation sort looks at sequences in ascending 
order of numerical value, sequences passing the test 
necessarily begin With three zeros. (The search is 
exhausted prior to reaching sequences which have a one as 
the first bit on the left.) then a Potential unique word 
(UN) sequence passes the autocorrelation test it is known 
that its ccrplement (CUH), transpose (TU) aná conplement 
transpose (CTU) will also pass the test. The latter were 
computed and stored along with the UW. For example, the 
following seguences of length 14 pass the autocorrelation 


test: 


Ug CUH TUR CTUS 
00000101100111 11111010011000 11100110100000 00011001011111 
00001011001111 11110100110000 11110011010000 00001100101111 


The cross-correlation condition is identical to that of 
autocorrelation on the head and tail cf the distribution. 
That is, the cross-correlation between sequences must be 
equal to zerc for the first and last three values. It is 
obvious from the example above that  cross-correlations 
between the UH/CTUW sequences and the CUN/TUK sequences do 


met meet this condition. This 1s true in general. Further, 
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because orf the symmetry inherent in the complementing 
Operation, (1) if UW1 cross-correlates well (matches) with 
CTUR2, then CUW1 matches with TUW2 and (2) if UW1 matches 
with UW2, then CTUW1 matches with CTUW2, etc. Therefore, in 
checking cross-correlation it is only necessary to check 
each UW against other UW's and CTUW's. Each set of unique 
words obtained by this process can be converted, by 
inspection, to a second set by taking the complement of each 
unigue word. This provides a dual flexibility in applying a 
set to a particular operating system with its own peculiar 


arrangement of fixed value bits preceding the unique word. 


CXUUNIQUE WORD SETS | 

AS explained above, the sequences fron the classes 
designated CUW and TUW need not be, and vere not, considered 
in the cress-correlation search for sets of unique words. 
Each set extracted from the UW/CTUW sequences has a 
companion within the CUW/TUW sequences, obtainable by 
inspecticn. eo uecurdcely -» befiect the 
conputational requirements, . the reported results which 
follow thus are drawn only from the former. That is, such 
quantities as “number of sequences passing the 
autocorrelation test!" and "number of sets of length X" refer 
to actual computer results. The numbers can be doubled by 
including the complemented sequences. Computer program 
descriptions are included in appendix D. 

Two results of primary interest were obtained for each 
value of N. First, the unigue word which cross-correlated 
well with the maximum number of sequences was reccrded, 
along with the corresponding matching sequences. These are 
referred to as "gated sets." These sequences may be used as 
unigue words in a gated acquisition system where it is 
necessary that only the reference uniyue word 
cross-correlate well with the other unique words. That is, 
as all detectors except the one set for the reference are 


gated, only the latter must contend with non-associated 
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unigue words during acquisition, and is thus the only one 
for which cross-correlation properties are of cencern. 
Second, sets of unique words were obtained which have the 
property that every word in the set cross-correlates well 
with every other word in the set. Sets of this type are 
necessary for ungated acquisition, and are referred to 
hereinafter simply as "sets." 

With h=2 and a fjve bit aperture width, no sequences 
passed the autocorrelation test for N=13 and below. At 
N=14, four sequences passed the autocorrelation test, but 
none of them could be paired in cross-correlation. At N=15, 
36 sequences passed the autocorrelation test. A gated set 
of length 13 and a total of 16 sets of length three were 
obtained. The total computer search time was 25 secs. 

At N=16, 150 sequences passed the autocorrelation test, 
and the largest gated set was of length 66. A total ci 56 
sets of length six were found. The computer search tise was 
3 min, 50 secs. 

The difference in computational requirements between 
N=15 and N=16 was the first indication of the magnitude of 
the time problem which would be encountered in dealing with 
higher values or N. At. N=17, the number of 
cross~correlation pairs increased to the point where it was 
also necessary to revise the final matching algorithm to 
prevent memory overflou, compounding the time problem. With 
444 sequences passing the autocorrelation test, a gated set 
of length 272 was obtained. Several sets of length nine 
were found using 4+ hours of computer time, at which point 
program execution was terminated. It is estimated that 
approximately 65 hours would be required for a complete 
search with the algorithm being used. 

It was obvious that it would be futile to continue the 
search tc higher values of N using the Same procedure. 
Attempts at developing a more efficient algorithm tere 
unsuccessful. The 32K on-line memory of the SDS-9300 was a 


serious ccnstraint, but even without it, the anticipated 
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execution time for the algorithms envisioned far exceeded 
the computer time that could be made available. A simple 
subterfuge was eventually employed to obtain partial results 
for higher values of N. 

The numter of sets of a given maximal length which were 
being found for higher values of N were far greater than 
needed in order to offer an acceptable degree of flexibility 
in selecting a particular set for use in an operational 
systen. Many of the sets differed only in the replacement 
of cne gr two unique words, and, because of the nature of 
the searching process, tended to be uncovered within the 
same time frame. When the program was terminated pricr to 
completicn cf a search, as at N=17, not only was it unknown 
whether or not sets of greater length than found actually 
existed, the ones that were found were for the mnmost part 
different ccrbinations of the same slightly larger set of 
unigue words. Given the fact that a partial search must 
suffice, it. seemed appropriate to effect a reduction in the 
number of seguences which passed the autocorrelation test to 
the point that the cross-correlation search could be run to 
completicr. Doing so again creates uncertainty as to 
whether cr not sets of maximal length have been found. 
However, given the excessive number of sets being obtained, 
there 1S some hope that some of them might still appear when 
the autoccrrelation basis is reduced. 

Increasing the effective detector aperture creates a 
more strict autocorrelation condition and thus reduces the 
numter of seguences passing the autocorrelation test. A 
nathematically artificial increase affects detector 
performance cnly in that the probability of false detection 
is reduced (unnecessarily) in the region immediately outside 
the real aperture. fhe effect of a varying aperture was 
first tested using N-15 and N-160, for which complete results 
had been previously obtained. For N-15, the maximal length 
set remained at three for aperture widths of three through 


AS he number of scch sets totaled 121 for a taree 





bit aperture versus six for a six bit aperture. For N=16, 
the maximal length remained at six for aperture widths of 
three through six bits, with corresponding totals of 724, 
345, 56 and 8. Four sets of length five were obtained for 
an aperture kidth of seven bits. 

Indeed, up to a point, the maxinal length is independent 
of aperture variation. Partial results were obtained for 
values cf N equalling 17, 18, 19 and 20 using this 
technique. The aperture was set high and then reduced in 
increments until search time became excesSive. The results 
so obtained again fell short of expectations in that much 
computer time was still required to complete the partial 
search. For a given number of sequences passing the 
autocorrelation test, it was expected that search time would 
be of the same order of magnitude for the different values 
of N. In fact, search time increased rapidly with N. For 
example, with a five bit aperture, 150 sequences passed the 
autocorrelation test at N=16 and the total search time was 3 
min, 50 sec. Wath a 16 bit aperture, 142 sequences passed 
HI EXutoccrrelation test at N=20 and the search was only 
approximately 1/30th complete when terminated after two 
hours. fhe imposition of a stricter autocorrelation 
criteria prcduces a basis of the same order as for a lower 
value of N, but the cross-correlation. properties of these 
sequences is such that many more combinations of match sets 
exist. ‘Zhat is, any given unique word cross-correlates well 
with many more sequences in the N=20 basis than in the N=15 
basis. 

Data on the numbers of sets found with h=2 for various 
values of N and various aperture widths is presented in 
Table 4-2; "4" is the aperture width and "NAUTO" is the 
numter of sequences passing the autocorrelation test. An 
asterisk denotes partial results obtained from an inconplete 
search. Table 4-3 contains similar data for h=t and bei, 


and serves to illustrate the effect of a varying threshold 
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on the numbers involved. A representative listing of 


maximal length sets ‚(for h=2) is contained in Appendix A. 





we vm, 





Gated Set Set Number of 
Maximal Maxinal Haxinal 

N R NAUTO Length Length Length Sets 
14 4 4 0 0 - 
15 7 1.10 2 2 D 

6 26 7 5 6 

5 36 13: 3 16 

4 60 21 3 T2 

3 78 22 3 121 
16 7 70 20 5 4 

6 110 47 6 8 

5 150 66 6 56 

4 190 84 6 345 

3 222 95 6 724 
17 9 56 21 6 2 

7 260 145 9% unkn 

5 LOD 272 9% unkn 
18 13 20 14 4 18 

12 46 26 6 63 

11 60 47 9 4 

10 142 102 10* unkn 

9 358 254 unkn unka 
19 15 36 26 5 89 

14 106 83 9x unkn 

13 204 157 unkn unkn 
20 16 142 113 10* unkn 

15 358 306 unkn unkn 


Table 4-2 
Sets Obtained For h=2 
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Gated Set Set Number of 


Maximal Maxımal Maximal 
h N.  NAJTO Length Length Length Sets 
1 10 7 4 2 2 5 
5 8 3 2 2 
11 9 8 2 0 = 
7 14 6 3 10 
5 30 13 4 20 
12 9 30 14 5 1€ 
7 DA 25 6 20 
5 18 qu 9 10 
13 9 92 54 10% unkn 
7 130 84 11* unkn 
5 182 120 unkn unkn 
3 19 5 0 0 0 = 
20 7 14 0 0 - 
5 16 17 y 26 





Table 4-3 | 
Sets Obtained For h=1 and h=3 


D. FALSE DEIECTION FROBABLLITIES 
The zrobability of miss-detection is, fcr constant h, 
the same for ali Sequences of a given length, 
Representative numbers are contained in Table 4-1. Given an 
acceptable ccmbination of length and detector threshold, it 
is the relative values of false detection probability that 
then determine whether or not a particular sequence, or set 
of seguences, 1S suitable for use aS a unique vord. The 
auto- and cross-correlation criteria were established with 
this in mind. As a point of final analysis, the equations 
of Appendix C were used to calculate the false detection 
prokapilities. Representative results are included in the 
foliowing paragraphs. 
1. Frobability of False Detection in the Random Region 
In an efficient TDMA system, the number of random 
(information) bits in a frame far exceeds the number of 
fixed bits. The majority of terms in the expression fcr the 


probability of one or more false detections per frame, 
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eguation C.10, are therefore of the form given by equation 
Gu, the probability of false detection when only random 
bats are in the detector. With h=2, 


eat, 


3.693 x 10 for N=15 


P {A (t}} 
r 


-— 


2.012 x 10 for N-20 


2. Crossccorrelation Bound 

When cne or more unique word bits, other than the 
one for which the detector is set, are contained in the 
detector, the cross-correlation criteria establish an upper 
bound on the probability of false detection (equations 3.4, 
3.6 and 3.7). These individual event probabilities can in 
turn be inserted into equation 3.8 to obtain an upper bound 
for the probability of one or more false detections (Ft) as 
any one unique word passes through the detector. Several 
such calculations were performed; for N=15 and N=20, the 


results were, 


= 3 
5.831 x 10 for N=15 


il 


Ft 
na X 


zu 
6.332 x 10 for N=20 


Thus F' is within an order of magnitude of the 
individual event probabilities in the random region. 
Because of the typically large number of random events, Ei 
in fact makes an insignificant contribution to the total F 
over the pericd of one frame. Given the fact that a 
cross-correlation false detection pulse is potentially 
periodic, an "insignificant" PI is precisely what is 
desired. 

3. FExatple System Calculations 

In crder to compute the probability of one or more 

false detecticns over a frame, it is necessary to have 


available the specific frame composition of the operating 
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system in guestion. A system was hypothesized for the 
purpose of obtaining a general feeling for the relative 
magnitude of the different factors involved. The number of 
Stations was set at three in order to simplify the 
calculations and permit selection of a set of unique words 
from those available at N=15. The fcllowing set was 
selected at random: 


(1) 000001100101111 
(2) 00010C100011111 
(3) 000101001101111 


This set possesses the following correlation 


properties: 
AUTCCORRBELATION CROSS-CORRELATION 
Bit) Rit) 

at 411 424) (3L t {21211} 432247) (3) 212) 
0 15 15 15 -15 0 0 0 
1 .9 9 7 -13 0 0 0 
2 7 / 7 ~12 0 0 0 
3 7 7 7 - 11 Q 0 1 
4 5 6 6 -10 0 1 2 
5 4 4 5 -9 2 2 2 
6 5 3 6 -8 y 3 3 
7 5 5 2 -7 5 4 2 
8 ü ü 3 -6 5 5 5 
9 2 2 3 =5 7 4 5 
10 1 1 2 -i 6 5 3 
41 0 1 1 -3 D 8 6 
12 0 0 0 -2 7 8 10 
13 0 0 0 -1 11 6 8 
14 Q 0 0 0 11 12 10 
+1 9 10 10 
+2 9 8 8 
+3 9 8 8 
+4 3 6 5 
+5 3 3 6 
+6 5 6 4 
+7 3 3 5 
+8 u 4 q 
+9 2 3 3 
+10 Z 2 1 
+11 1 1 1 
+12 0 0 0 
+13 0 0 0 
+14 0 0 0 


The hypothetical system is patterned after the 
operational system described in Ref. 3. The transmission 


rate is 60 Mbps (4-phase), the frame longt is 750"usec2lmbhe 
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guard band is 24 bjts in width and there are 60 fixed bits 
in the preamble for carrier and bit timing recovery. The 
latter are arbitrarily assumed to consist of a sequence of 
30 i's for carrier recovery, followed by a sequence cf 30 
alternating O's and 1's for bit timing recovery. Bits are 
clocked in from the left, with the detectors set tc the 
unigue wcrds as listed above. 

It is first necessary to determine the nunker of 
random versus fixed detection events per frame. The product 
of the transtission pit rate and frame length gives the 
number of bits in one frame, n=45000. The beginning of each 
station's burst is composed of 75 fixed value bits (recovery 
plus unique word). When one or more of these bits is in the 
detector, the detection event is, by definition, "fixed." 
Therefore, for each station's burst, there are 75 fixed 
detection events as each bit, from the first to the last, is 
clccked intc the detector; plus an additional i5 fixed 
events as the last bit, and those immediately preceding it, 
are clocked out of the detector. The total number of fixed 
detection events per frame is thus equal to 270 (3x90). The 
number of random events 1s 44730 (45000~270). 

With this number of random events, the upper bound 
on the probability of one or more false detections per frame 


in the random region alone is (equation 3.10), 
r3 
F< rP{A(t)} = 44730(3.693 x 10 ) 
r 


F s 165 


ech is net a particularly enlightening result for a 
Miptobability upper bound. It is indicative, however, of the 
fact that a false detection is almost certain during 


acquisition. By way of comparison, the Ref. 11 expression 


qu 











for the probability of no false detection in the random 
region gives, 


44730 
1- F = [1 - P{A(t)} ] = 1.3 x 10 
I 


Therefore, it can be concluded that many false detections 
will cccur during the acquisition phase, and that, fcr the 
System design hypothesized, little can be done to improve 
the situation, snort of using some type of acquisition gete. 
This is generally true for the typical prototype TDMA 
systems which have been developed. No attempt was made to 
compute the expected acquisition time, a subject treated 
extensively in Ref. 11. 

If an acquisjtion gate is employed, the detectcr is 
active fcr a maximum of 2B=48 detection cycles pricr to the 
receipt of the associated unigue word. All of these are 
tixed detection events. At each cycle, the autocorrelation 
values can be added to the number of matches occurring in 
the recovery bits to obtain the number of total fixed 
matches, k.e Equation C.8e can then be used to o taaan 


upper bound cn Fs. Doing so, for detector number (2), 
-11 
ps Ego x 
| 10 
|o 02503556 09 frames = 217 days 
For detector number (3), 
- y 
ESA ISO 
F > 1669 frames = 1.25 sec 
The mean time between false detections (fF) is more 
than acceptable for both, indicating a high probability of 


achieving acguisition in just one pass of the unigue word. 


Acquisition cf the reference unique word, assumed to De 
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number (1) here, must still be effected without benefit of 
gating, and thus remains a potentially time consuming 
process. 


LO 





V. CONCLUSIONS 


This research effort was undertaken for the express 
purpose of finding unique word sequences which would be 
suitable for use in a TDMA satellite communication system. 
The correlaticn properties of potential unique words vere of 
primary ccncern in this effort, criteria being established 
which wculd insure a low probability of false detection. 
Numerous sets of unique words were found using these 
criteria. 

The results obtained can not, unfortunately, be applied 
to any given operational TDMA system without some additional 
consideraticn of false detection probabilities. In an 
attempt to generalize findings, it vas necessary to ignore 
the fact that each system has its own peculiar set of fixed 
value bits, in addition to the unique words, which are 
transmitted iu each siaition's burst. These fixed value bits 
are in fact an extension of the unique word sequence and 
produce a cross-correlation distribution of their otn when 
passing through a detector. The cross-correlation may very 
well be such that a set of unique words selected at random 
from this work would prove unacceptable. Fortunately, the 
number of sets available of any given length, together with 
the existence of the complementary opticns, provides a 
flexibility sufficient to cover any foreseeable situation. 

The length of the sets found fell short of expectations. 
To begin with, in setting cross-correlation criteria based 
on worst-case possibilities for the number of detector 
agreements, unique words which might actually be useable in 
a particular system were excluded. That iS, the number of 
agreements produced when non-unique word fixed bits are 
stepping through the detector can seldom be expected to 
approach the worst-case limit. This fact was noted at the 


Butset. however, and accepted in order to generalize the 
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results as much as possible. The real disappointing factor 
was that, even with the strict correlation criteria, a 
computer search could not be designed tc uncover with 
certainty the sets of maximal length for values of N in 
excess of sixteen. 

There is a need for an efficient algorithm for comkining 
unigue wcrds into sets possessing good cross-correlation 
properties. As stated in Ref. 11, the problem is a puzzling 
one. It is a problem which can undoubtedly be solved only 
by adopting some type of short-cut technique, wherein but a 
few of all the multitude of different patching combinations 
are considered. The ability to predict the maximal length 
for a given N would be of tremendous benefit in this regard. 

The ultimate achievement would be the discovery of an 
underlying structure which enabled one to extract unique 
word sequences without the need to resort to brute force 


computer search techniques. 
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SETS 


Selected sets of unique words, representative cf the 


N=15 is 
unigue words are listed 
valid 


many found, are presented herein. The listing for 


complete. In the interest of space, 
All 


threshold setting of two. 


in octal fcrmat. sets are for a detector 


A "gated set" is one in which the first sequence listed 


(the supposed reference) cross-correlates well with all the 
sequences which follow. A "complete set" is one in which 
each Sequence cross-correlates well with every other 
sequence. 
RIO 
MAXINAL LENGTH GATED SET 
01457 01347 02067 02327 02367 02467 03067 
04437 05157 04157 04657 04517 04117 
COMPLETE SETS 
un 00627 02067 02217 (IN - 01667 02367 02717 
(2) 06627 02217 02327 (10) 01667 02367 02637 
(3) 06627 02217 02367 (11) 01667 03367 02637 
(4) 01317 01667 02327 (12) 05477 04757 03557 
(5) 01357 04437 05157 da zur, 5903557 (05 159 
(6) 01457 04437 04157 (14) 05477 03557 04157 
(7) 01457 04437 04117 (15) 03137 04437 05157 
(8) 01667 02327 02637 (16) 04437 04157 03057 
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MAXIMAL LENGTH GATED SET 


002637 


OO SO COOoODd0O 
web chte EE EEN 


OOO2EFDEEFELIN> 


COMPLETE SETS 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 
(13) 
(14) 


061457 
001667 
002367 
002467 
002637 
002637 
602637 
002717 
003057 
003267 
003557 
004567 
004417 
005317 


ZSANINWEN SUN 
ANNIE NND OD 


Ja SS as ay Sy 
On EE EE EEN 
ma mm wim am Lë Le Leg Lo La 
O a MNO m O O E EUN Iam 
3 VO WN UU o Y EM 
æ YN DN Oy LY Da NE DES 


002317 
002637 
005.167 
002637 
003267 
003267 
004657 
005317 
004567 
004367 
011537 
005067 
003137 
011077 


N=1€ 


- 


SNR SUR Sd SS 
E EE EE EE 
Ome OOOO OO 
eNO & » 06 ElyN a 
ONO ON Uy Oy Y MON 
eet ECKE EE EE 


005317 
004657 
0006467 
003267 
004367 
005167 
006137 
010337 
005067 
004717 
005637 
006267 
004437 
010337 
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VM SOS RURSUS SU 


011237 
005 167 
005477 
004367 
006467 
006467 
011237 
003137 
006267 
005467 
010557 
006137 
005157 
003137 


SSOS999595 2595 95 
—OOOO-oOocoooco 
DANIAN Mme DH $ GI Woh 
NR) a Ss SJ ONO a LA 
wm LIDL) LA LA UN wet CN POND 


010657 
006467 
006137 
006137 
006137 
006137 
010357 
010657 
004437 
003137 
011657 
004437 
010157 
010657 


aju Se Sal Se Se Sal Se Se J «— 
EE 
OO 2.» OOOOO 
RDNS ENDS de GON) 
Ge e] o Ly NO O ao NY GO 
CANN CI LU CRIN A 


011517 
006157 
006257 
006257 
006257 
006257 
011457 
011517 
005157 
005157 
011317 
005157 
011517 
011517 


YYYY 








cu 
= 


La 


CUNELETE SETS 


012467 


007327 


011767 


024437 


007267 


007157 


007567 


sus 


COMPLETE 


004767 
021667 


Sd 





SereLETE SETS 


SES TER SETS 
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DETERMINATION 


In the absence of transmitter burst position corrections 
and under zero doppler conditions, a given unigue word will 
be received every T seconds. (Other factors, such as clock 
stability and logic jitter, affect this periodicity; 
however, they are jnconsequential compared to potential 
doppler effects and will be ignored.) Burst drift is 
defined as the doppler induced, real-time shift of a burst 
from its periodic position. 

Burst drift rate is a function of the relative velocity 
between transmitting station, satellite and receiving 
station. Interest here is in determination of the maximum 
drift rate that can be expected. First, let us consider the 
motion of the typical communication satellite in an 
equatorial, gecsynchronous or 
remains at the same altitude over a given point on the 
equator. In reality, it is physically impossible to achieve 
and maintain a perfectly synchronous orbit [ 12]. Reference 
13 derives equations which can be used to estimate the 
worst-case relative velocity between a geosynchronous 
satellite and a given point on the earth's surface. The 
orbit eccentricity and inclination are the two primary 
factors which influence the satellite's moticn. 

A non-circular orbit results in a sinuscidal variation 
of the satellite altitude. To a first order approximation, 
the resulting doppler shift is independent of earth station 
coordinates. If eis the eccentricity of the orbit, the 


maximum relative velocity is approximately, 


v = 12100e Kn/hr = 3361le m/sec 
e 





Using .0004 from Table 1 of Ref. 14 as a typical 

eccentricity for later generaticn communication satellites, 
v = 1.3 n/sec 
€ 

When the orbit is inclined from the equatorial lane, 
the satellite moves in a figure-8 pattern, elongated north 
and south. The north-south motion  predominates and is 
Sinusoidal in nature. If i is the orbit inclination in 
radians, the resulting maximum relative velocity as seen by 
a Station at latitude L can be approximated by, 

v. = 2092isin(L) Ku/hr = S58lisin(L) m/sec 
i 

Satellite visibility extends to approximately 80° 
latitude. At this Jatitude, the satellite appears just 
above the horizon and relative velocity due to orbit 
inclination is a maximum. (it is in fact highly unlikely 
that a staticn would attempt to use a satellite below 10° in 
elevation.) Thus, assuming a maximum inclination cf 3° 
(-052 rad) and L=80°, 

v = 29.7 m/sec 
1 

A pessimistic estimate of the maximum relative velocity 
due to satellite motion alone can be obtained by comktining 
the effects of orbit eccentricity and inclination, giving, 

v =v + v = 31 m/sec 
S e i 

Next, the contribution of earth staticn motion must be 
considered. Reference 13 analyzes the headway, roll and 
pitch components of ship motion and concludes that the 
maximum pcssible relative velocity contribution is, 


v = 38 m/sec 
g 





For a shipboard platform,then, the maximum relative 

velocity between ground station and satellite is, 
V=v +v = 69 m/sec 
g S 

With representative numbers for satellite and surface 
platform motion in hand, it is next necessary to consider 
how they right combine to produce a shift in the position of 
a received unique word. The case of a station receiving its 
own unique word (transmitted up to the satellite and back) 
is a simple two-statjon relative motion problem and the 


weli-kncwn dcppler equation applies. 
d = 2vf/c (B.5) 


where f is the transmitted bit rate and c the speed of 
light. 

The case of a station receiving another's unigue werd is 
a three-station motion problem. Maximum doppler occurs when 
the maximum reiative velocity exists between both ground 
stations and the satellite simultaneously; that is, when 
maximum  up-link doppler is combined with maximum down-link 
doppler. Such a situation is produced when both ground 
Stations are at the same 809 latitude and both are either 
opening or closing the satellite at the maximum rate. For 
this case also, then, with maximum relative velocity, v, 
existing between both stations and the satellite, equation 


B.5 applies. 








DETECTION PROBABILITIES 


Expressions for the probability of missed and false 
detections are derived in the following paragraphs. These 
expressicns depend somewhat on general detector operation, 
but not to the extent that their application is restricted 


narrowly to any one particular systen. 


Roa nLoS- DETECTION PROBABILITY 

If the bit error rate (BER), i.e. the probability that a 
received bit will be in error, is p and the number of bits 
in a unique word is N, then, assuming errors occur 


independently, the probability of a unique word being 
l 
received error-free is (l-p) . With the detector threshold 


set at h, a unique word will be detected if it contains h or 
fewer errors. The probability of detection, P, is thus 
obtained by summing terms in the binomial distribution which 
correspond tc the probability cf no errors, the probatility 
of one error, up to the probability of h errors. Using 
C(N,i) to represent the binomial coefficient N!/i:? (N-i), P 
is given ty, 


h n i 
E SC C(N,1)p (1-p) (C.1) 


The miss-detection probability (the probability of 


receiving a unique word with more than h errors) is given 


by, 
N Ee SE 
Q = 1-P = 2 C{N,i)p (1-p) (C.2) 
i=h+1 
The mean time between miss-detections is equal to the 
reciprocal of the expected value of "number of misses per 


wat time," As there is only one possible miss-detection 
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event per frame, the expected value is equal tc the 
miss-detection probability, Q, and the mean time between 
miss-detecticns is, 


-1 -1 
R= Q frames = TQ secs (C.3) 


B. PROBAFILITY OF FALSE DETECTION 
The following notation will be used in developing 
expressicns for false detection probabilities: 
‘A (1) “RK false detection event at time t. 


Used interchangeably with "A ." 
t 


E (A (t)] The probability of A(t). 
E (t) No false detection at time t. 
Used interchangeably with "B ." 
t 


E (B(t)] Ihe probability of Bt) . 
t An integral time variable, the basic unit 


of which is one bit cycle (t=ib). 


On any given bit cycle, other than the one in which the 
actual unique word is completely contained in the detector, 
a false detection pulse will be produced if N-h or mere of 
the N bits in an active detector match with the detector 
setting (R2h). The probability of this occurring depends on 
whether cr net the value of transmitted bits are fixed or 
variable. The beginning portion of each burst is composed 
of fixed-value bits (timing recovery and unique word), while 
the remainder of the burst contains bits of variable value. 
The probability of one of the latter being transmitted as a 
(01 or a *'1* depends On the type of informaticn transmitted. 
As a first order approximation, it will be assumed that the 
variable bits are in fact random, taking on the value ‘0! 


and '1* with equal probability. With random bits alone in 


è 
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the detector, the probability of a false detection cn any 
given bit cycle is, 


h . i N 
P(A(t)). 7,2, CQ, i) (1/72) (172) 


N h 
= (1/2) = C (N, 1i) (C.4) 


The probability of false detection is not so easily 
computed when there are fixed value bits in the detector. 
"Fixed value" is used to denote those bits in the preanble 
which are transmitted with the sane value ('*0' or '1') in 
each frame. Because of the presence of transmission errors, 
these bits may actually be received with the opposite value, 
the probability of which is p. An expression fcr the 
probability of false detection can be derived based on the 
number of fixed bits in the detector and the number that 
matched with the detector setting when transmitted. Such an 
EUEeSsicH uust accoHut ror ali tie possible combinaticnes ef 
errors within the fixed bits and as such will include 
additive terms containing higher powers of p. With p<<1, a 
good approximation can be obtained by discarding all terms 
which contain other than the lowest power of p. Therefore, 
rather than derive a complex general expression from which 
an approximation can be drawn, the condition p<<i will be 
invoked at the outset in order to simplify the derivation. 

Suppcse there are m fixed bits in the detector and that 
k of them matched with the detector setting when 
transmitted. A false detection will occur when the number 
of agreements within the fixed bits, plus the number of 
agreements within the remaining N-m random bits is greater 
than or equal to N-h. With p««1, interest can be focused on 
Muse particular combinations which require the fewest 
number cf transmission errors in order to produce a false 


detection. Three separate cases will be considered. 





CASE 1. The number of fixed agreements, k, is 
greater than or equal to N-h. A false detection is assured 
if them fixed hits are received error free, and the 


probability cf false detection can be approximated by, 
n 
P{A(t)} SE P EE (C. &a) 


CASE 2. The number of fixed agreements, k, is less 
than N-h, but enough random bits (N-m of them) remain to 


produce a totai of N-L or uore agreements. That is, 
kit aN eh 
Or, 


k + h 


1V 


m 


A false detection will be produced when there are no errors 
in the fixed bits and N-k-n or more agreements withir the 


random bits; or, equivalently, when there are 
Se eh ea n 


or fewer disagreements within the random bits. The 
probability of false detection can thus be approximated by 
multiplying the probability of receiving the n fixed bits 
error free by the probability that the N-m random bits will 


contain kth< or fewer disagreements. 


P {A(t = 
at, 


En k+h- 


m N-n | 
(i-p) (1/2) 2 Com) m<k+h<N (Ca Eb) 


— 


CASE 3. The number of fixed agreements, k, plus the 
number cf random bits, N~m, is less than N-h; i.e., kth<m. 
A false detection can not occur unless there are errors in 
the fixed tits. The fewest number of errors are reguired 


when all of the random bits and all of the original k bits 
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match with the detector setting. In this case, enough 
errcrs must cccur in the remaining m-k fixed bits to produce 
N-h agreements (the minimum required for false detection). 


The number of errors required is, 


a A (k+ N- n) 


=m- k - h 


The protatility of false detection can thus be approximated 
by taking the product of (1) the probability that the N-a 
random tits will all match with the detector setting, (2) 
the probability of receiving the original k matching bits 
error free and (3) the probability that the m-k fixed bits 


will contain 'a' errors. ` 


E. 


m k a n-k-a 
en., = (1/2) (lap) C(m-k,a)p (1-p) 


N-m a m 
(1/2) C(m-k,a)p (1-p) 


D n-xcd k*h 
C(m-k,m-k-h) p (1-p) ; k+h<u 


(C. Ec) 


N 
(1/2) 


If all the bits in the detector are fixed, the 


probability cf false detection in any given bit cycle is 


given by equations C.8a and C.8c with m=N. 


N 
Auer zielen); k+h2N (Ce €d) 


kth 


-k-h 
(1-p) S k+th<N (C. Bei 


N 
P{A(t)} - C(N-k,N-k-h) p 


—— ee ee diu wee ar. ae age vm: ffr cu caue m 


Fouations Cp and C.8 give the probability of false 
detection for a single detection event within the frame. 


Calculation of the probability of false detection over 
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several consecutive bit periods, ‘specifically over one 
frame, is complicated by the fact that successive detection 


events are nct independent; i.e., 
P(A(t),AÀ(t* 1)]) *P(A(t)]) P(A(t*1)). 


That this is true can be seen by noting that successive 
detection events are generated by shifting the contents of 
the shift register by one bit. If a false detection does 
(or does nct) occur at one step, possibilities for the 
"contents ol the shift register are necessarily restricted to 
certain combinations of bit values. These bit values carry 
over to the next detection event and reduce the sample space 
accordingly. As a simple, extreme example, consider the 
case of the unique word 101010. The probability of false 
detection with random bits in the detector and a zero 
threshold setting 1s, from equation C.4, (0.5) . However, 
given that a false detection occurs {the detector contains 
101010), the detector will contain X10101 on the next step 
and the probability of false detection is zero. 

Dependent joint probabilities are best dealt with by 
locking for, and finding, grounds for the assumpticn of 
independence. Schrempp and Sekimoto [11] based their 
independence assumption on the previous work of Hill [14]. 
The validity of this assumption is nearly impossitle to 
access without making lengthy comparative calculations for 
unique words of different length, a task of monurental 
proportions. Rather than attempt to justify such an 
assumpticn here, joint probabilities will be expressed in 
terms of an upper bound. 

There are n detection events per frame during 
acquisition, n-1 of which may produce false detections. 
Defining t-0 at the point when the associated unigue werd is 


completely contained in the detector, the probability cf no 
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false detection during a frame period is given by, 


E 
1 


P{E ,B ,+..,B 3 
{ 172" d -— (3:9) 


PÍE JP(D /B )JP(B /B B ) eeoeP(B B b eeeBp 
d : ai d y 1 2 i i 2 En 


Cr, from a different viewpoint, if F is the 
probability cf one cr more false detections per frame, 


H = 1 - P 
SEN n-1 n-1 
F= Ż Pea SE PO eh 
1-1 1 i-1 "Ei 1293 
n-1 n-1 n-1 
+ 


161 Zi 5j 


t P(A ,A ¿90e,A 
( es i m 


For straight-foruard calculations, F is obviously 
orders of magnitude more difficult to obtain than is H. 
However, its form does offer possibilities for reducing the 
complexity cf „the problem. Although successive events are 


not independent, events separated by nore than N detection 


Cycles are - - after N detection cycles, all bits frcm the 
First event have been shifted out of the detector. For 
HON, 


EE 
i j i j 


and similarly for higher order joint probabilities. If 
P{A(i)}<<1, these product terms are Small compared to the 
individual terms from the first summation in F. However, 
mere the dependent joint probabilities it is quite possible 


that 


P{A ,A } = P{A JP{A /A 3 
5 i JP 
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is of the same order as PIA }. The same can be said for 
i 
higher crder terms, but with a decreasing expectation cf its 


actually  prcving true. It is therefore not possikle to 
discard as insignificant all of the terns involving joint 
probabilities. The most that can be said is that, 


P(A )2P(A ,A J2P(A ,À ,A ) e o e 
i i 227 i j K 


And that, 


n- 1 
Ben PÍA ] 
1-1 1 


The false detection probabilities Pf{A(i)} are 
obtainable from equations C.4 and C.8 as appropriate. For 
the randcn detection events (eguation C.4), 
P{A(1)}=EFE{A(j)}: for the fixed detection events, tae 


prokabilities must be computed individually. If there are r 


random faise detection events and  f-n-i-r fixed false 
detection events per frame, 
n~- 
F = 2 P(A } 
max 1-1 1 
= EP{A(t)} + 2P1A(t)) (C. 10) 
r £ f 


When the detector is gated during acquisition, as 
described in Section II.C, the detector is active for only 
2B+1 bit cycles about the nominal position of the associated 
unique  wucrd. The true detection pulse may occur anywhere 
within the aperture, depending upon the amount of 
displacement of the unigue word from its ncminal position, 
and the maximum number of false detection events is equal to 
E The probability of one or more false detections in the 
aperture prior tc the receipt of the unique word is 
therefore of the same form as equation C.10, with the number 


of events equal to 2B vice n-1. 
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In the steady-state mode, the aperture width is 2D+1 
and the maxinum number of false detection events is equal to 
2D. 

In Summary, equation C.10 provides a bound on the 
probability of one or more false detections over one frame 
period, or, with a change of limits, over an acguisition or 
steady-state detection aperture. This bound is useful in 
estimating the freguency of false detections and in 
assessing the relative merits of different unique word 
comkinaticns. 

AS an aSide, it is interesting to compare the randon 
event term of equation C.10 with the results of Ref. H, 
obtained by assuming independent events. For randcm bit 


detection events, it was concluded that, 


E O SS P{A(t)} ] 


Expanding, 
2 
F=1- [1 cP{A(t)} + C(r,2) P{A(t)} - eee] 


For E{A(t)}<<1, a basic assumption cf Ref. 14, an 


approximate expression for F is, 
EOD AD 


Which is identical to the randcm bit term in equation  C.10. 
2. Mean Tine Betzeen False Detections 

A false detection, both in the acquisition and 
steady-state modes, deactivates the detector and sets an 
aperture for the next frame. Thee probability of cycling 
through several successive detection events without any 
false detections is therefore of primary ccncern. Fron the 
opposite viewpoint, it is the probability of one or more 
false detecticns over a given time span which is meaningful, 
not the probability of false detection at any given step. 


This distinction is important in deciding upon what tc use 
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aS an appropriate measure of “mean time between false 
detections" (MTBFD), and in turn interpreting its meaning. 
If the detector were not deactivated by a faise 
detection pulse, several false detections might occur over 
the pericd of one frame. The MTEFD would then be equal to 
the reciprocal of the expected value for "number of false 
detections per frame." With detector deactivation, only one 
false detection per frame can actually occur and the MIBFD 
is affected accordingly. In this case, for computational 
purposes, the probability of a real, physical false 
detection pulse occurring somewhere within the frame is 
equal to F, the probability of one or more false detections 
per frame. The MTBFD is therefore equal to the reciprocal 


of F and can be bounded by substitution from equation  C.10. 


" -1 -1 
& = F > F frames (C. 11) 
max 


65 











The ccmputer program was written, and dekugged, in three 
Separate functional segments and then patched together to 
form one continuous program. These segments will be 
described and listed separately. The last segment listed is 
the basic cross-correlation matching routine. Several 
different revisions were employed in an attempt to overcome 


time and memcry limitations. 


Arrays 


GUW (ICIN) —- Sequences uhich pass the autocorrelation 


test. 
GCTUBWB(IDIM) -~ Complement transpose of GUW. 
GTUR (IDIN) - Transpose of GUM. 
GCUN (1DIM) - Conplenent of GUW. 


MTEME(IDIN2) - Used for temporary storage of good 
crces~correlation matches. 


INCR ( ) - Used to store pointers needed in the 
crcss-correlation E process. Dimension greater 
than or equal to maxımal length set expected. 


BRICH(IDINT,IDIN2) - Hatrix in which good 
Cross-correlation matches are recorded. 


HCROSS(21DIM) —- Sequences from autocorrelation sort 
which are checked for cross-correiation (the GUH and 
GCTUN subsets). 


Equivalence Statement - The first half of MCROSS is 
Stcred in same memory location as GUW( ); the second 
half in same location as GCTUW( ). 


66 





Variables Common to All Segments 


UW - Sequence under test. 


SUH - Shifted version of UW. A sequence is shifted cne 
bit at a time and, at each step, added mod-2 to UW. 


TUN - Transpose of UW. 

CTUW - Complement transpose of UW. 

CUH - Couplement of UH. 

N — Length of sequences being tested (Input). 
IEPS - Threshold; h of thesis main body (Input). 


IALFÀAÀ ~ Aperture parameter; equals 2D (Input). 
NUW RE of good UW's obtained in autocorrelation 
Sort. 


NAUTO ~ Number of UW's and CTUW's, combined, which are 
checked for cross-correlation; equals 2NUW. 
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arjables 


—— -— 800 --—— 


K - 5mpallest integer equal to or greater than N/2. 


LAST,LASI1 - Intermediate skip points in the sequential 
checking process. Sequences which must be checked 
tend to occur in blocked subsets, with two redundant 
Sequences separating each. 


I21 - Next to last sequence which must be checked. 


122 > Last sequence which must be checked to insure all 
possible non-redundant sequences have been considered. 


Mere Rinary constant consisting of N 1's, Eight 
justified. Used to obtain the couplement of an N DIE 
seguence. 


IFLAG - Set to one vhen last sequence being checked. 
Keys the termination of autocorrelation search. 


JFLAG - Set to one when autocorrelation condition 


requires SR eo Autocorrelation computation for a 
particular UW is complete when JFLAG=IEPS. 


-* 


I - Shift variable. Equivalent to time variable of 
thesis main body. 


ICOUNT - The autocorrelation value at step I. 


The sequences which pass the autocorrelation test, and 
the corresponding CUN's, TUW's and CTUW'S, are output in 


side-by-side decimal and octal format. 
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- Number of sequences with which sequence being 
tested has cross-correlated acceptabiy. Used for 
filling NATCH matrix. 


M2HIGH - Highest nunber of sequences with which any one 
Sequence haS cross-correlated acceptably. Used for 
MATCH matrix output. (This is not the maximum number 
of cross-correlation matches existing. Correlation is 
Checked, and matcnes stored bur once, commencing with 
the first sequence in,the AUTO-C set and working up. 
Thus H2HIGH does not include matches between the 
Sequence being tested and previous sequences in the 
AUTO-C set.) 


Other variables are defined under "Arrays" and 
Pconacn. " 


(1) The UW and CTUW sequences from the AUTO-C set which 
are being cross-correlated are oat oit in decimal fcrmat, 
together with an index number which denotes the order in 
Which the sequences are tested. 


(2) Tbe matching matrix (MATCH) as output in 
two-dimensional form. The elements of MATCH are the 
index numbers SEE Eon to the .AUTO-C “get of 
sequences. There is one row for each AUTO-C Uk (See 


Section IV.B). The first element cf a row is the index 
number of the UH, and the remainder of the row contains 
(higher) index numbers of the Uz's and CTUW'S with which 


it matches. 
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C HAVE WE CHECKED UW AGAINST ALL HIGHER ORDER SEQUENCES? 
IF{M.LE.NAUTO) GO TO 110 

C FINISHzD WITH UW. STORE MATCHES. 

120 M2=N2- 


esch 


DO 125 K=1,M2 
MATCH Lyk) =HIEMP (K) 
125 MTEHP (K)= 
C STORE NUMBER OF UW MATCHES IN FIRST UNUSED ROW ELEBENT 


MATCH (L,IDIM2) =M2 
C IS M2 HIGHER THAN FOR ANY PREVIOUS UN? 
IE (2.61. N2uIGH) M2HIGH=M2 


TIME TC GET NEW UW. 
L=1+1 


C 

C ABE WE FINISHED? 
IF (L.LE.NUR) GO TO 105 

E SORT COMPLETE. OUTPUT MATCH. 


met 


TE (6,655) 
137 I 
TE ( 


DO E 
137 WRITE(6 666) (MATCH (I, J) ,J=1,M2HIGH) 
S GO 360 
C GC TO HATCHING ROUTINE 
C DCES PAIR FASS ZERO TEST AT STEP I? 
140 | TF (ICOUNT.GT.0) GO TO 117 
JFLAG=dJFLAG+ 1 
C ARE WE FINISHED WITH RIGHT-SHIFT? 
IF (FLAG. .LT. IEPS) GO TO 115 
JFLAG=0 
SUW-MCAOSS (1) 
C START LOCE FOR COMPUTING LEFT-SHIFTED CROSS-CORRELATION. 
C THE IN-LINE SSGHENZ INITIALIZES N-I AND I COUNTERS, LOADS 
C SUw INTC A~REGISTER, SHIFTS SUN LEFT ONE BIT, STORES THE 
C RESULT FCR LATER USE, COMPUTES MOD-2 SUM OF UW AND SUH, 
C SHIFTS THE MOD-2 SUM RIGHT I BITS SO THAT THE N BITS OF 
C INTEREST ARE RIGHT-JUSTIFIED IN THE A-REGISTER, COUNTS 
C NUMBER OF ZEROS IN FIRST N-I POSITIONS (FROM RIGHT) CF 
C THE MOD-2 SUM. 
145 I-1-1 
ICCUNT=0 
S LDA N 
S SUE I 
S ADD 077677177 
S RCH 070100,1 
S LDA i 
S ADD 2077617111 
S RCH 070100,2 
S LDA SUN 
S SHIFT (25001 
S STA SUE 
S ECR UW 
S SHIFT 01001 
S BhX $-1,2 
S SHIFT 01001 
S SKB 20400000000 
S MPC ICCUNT 
S BRX $-3,1 
IF(I.GE.N-IEPS) GO TO 150 
IF tICOUNT.LT.N-IEPS-I) GO TO 145 
eo so 177 
EU IP(ICOUNT.GT.0) GO TO 117 
JFLAG-JELAG* | 
IF(JFLAG.LT.IEPS) GO TO 145 
C HAVE A COOL HATCH. STORE IN MTEMP 
MTEMP (£2) = 
H2=M2+1 
C DIMENSION CVERFLOW CHECK 
1E (2-17. Er 12-1) GO T9 117 
HRÍTE (6,665) L 
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RITE (6, 670) 
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ariables 


HOST - Length of largest set found at any given point in 
Search. Smaller sets, when found, are not output. 


LTEST - Sets the last row element which must be looked 
at Lased on current value of MOST. 


MATI ~ An element from the matrix row which has MAT2 as 
the leading element. 


MAT2 - The last element of the current set being formed. 


MAT3 - An element from the matrix row which has one of 
the set members as a _ leading element. Once it is 
kncun that a potential addition matches with the first 
andoJast elements of the set, HAT3 is Used to Check if 
it matches with the other elements as well. 


Other variables are defined under "Arrays" and "Common." 


Match sets are output, one line at a time in decimal 
format with the autocorrelation index preceding each element 
of the set. 
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Progras 
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C HAVE NEW ALDITION TO SET 

GO 10 355 
352 MATZ=MAT2-NUNW 

DO 354 MM1=MM,M1 

LL=2 
353 EE L 

TE IER OR.NAT3.GT BAT) GO TO 325 

IF (MAT3.EY.MAT2) GO TO 354 

LL=LL+1 

GOTTO 353 
354 CONTINUE 

MAT2=MAT2+NUW 
C STORE NEW ADDITION AND SET POINTER 
359 M=M+1 

HTEMP (M) -MAT2 

INCH { ed 
C START LCOKING FOR SEQUENCES WHICH MATCH WITH FIRST AND 
C (NER) n ELEMENT OF SET. 

R= 

IF (MATCH (I[,K).GT.NUW) ,LTESI=MATCH (MATCH (I,K)-NUW, 

1IDIE2) +M+1-NOST; GO TO 320 

LTESI=BATCH (MATCH (1,K) ,1D1M2) +M+1-MOST 

GO TO 320 
OR TORE FIRST THO ELEMENTS OF 2-ELEMENT SET AND SET 
C PCINTERS. THEN GO TO POINT WHERE NEW ELEMENT IS ADDED 
TO SET 
365 NTEMP(1)= 

NTENP {2} -MATCH (I,K) 

INCK(1)= 

TNER(2)=K 
2 GO TO 355 
C HAVE EXHAUSTED POSSIBILITIES FCR ADDITIONS TO CURRENT SET 
O SET LENGTH OF INTEREST? IF SO, OUTPUT; IN ANY CASE, 
Beh UE TO NECAD- TO-LAST ELEUENT CF Sen DE ee oe EPOR 
C FURTHER POTENTIAL ADDITIONS. 
370 IF(M.GE.MOST) GO TO 375 
372 IF(M.EÇ.2) GO TO 405 

GO TO 380 
375 MOSI=M 

WRITE (6, 205) (HTEKBUNT) „BEROSS (NTEAB (MTY) MI 1,MOST) 
380 ee Nek (1-1) 

J-INCH (M) +4 

IF (MATCH (I, H GT.NUW) GO TO 450 

HAT ISHATCH (I, KY 

MAT2-MATCH(I,INCAB(M)) 
385  L-L*1 S 

IF (NATC H(MAT1,L).LT.MAT2) GO TO 385 

M-H- 

ERRESTD-MATCH(MATCH(I,K),IDIM2)THMtT1-HMOST 

GO TO 327 
405 K=K+1 

IF (NATCH (I,K LES «GT. BE ,LTES T= MATCH (MATCH (I,K) -NUR, 

ED. H2) ums e 15 

LIEST= HATCH (HATCH (Lo E IDAZ tiri MOST 

GO LO 
C REMAINING STATEMENTS HANDLE CHECKS WHEN MATRIX ELEMENT IN 
C QUESTION IS FROM CTUW AUTO-C SUBSET 
415 MAT3=MATCH (1I,K) -NUN 
420 L=1+1 

TE (L +61. LTEST) CORTO T 

MATI=MAICH (HAT3,L) 

TF (HAT 1, EO. 0) GO TO 370 
425 MAÍ2=MATCH (1,3) - NUM 

IP(MATZ.LE.0) GO TO 370 

TF (BATLLE? MAT1) GO TO 430 

GO TO 4 
430 IF{NAT1.50.MAT2) „MAT2=MAT2+NUN:;GO TO 340 

Usa 

CGO 10 425 
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450 
455 


JO 


Ee ,K) -NUW 

BAIZ-HATCU I,INCR(M))-NUN 

=, + 

IF (HATCH (HATT, TI STI ADOOS 
LTEST=MATCH (MAT1,IDIN2) +M+1-MOST 

GO 10 415 

ES 590p CROSS-CORRELATICN SETS'//) 
FORMAT 101,5 (1X,7 (I4,!) ©, 1X,16,5X) /) ) 

E 
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